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ABSTRACT 

Context. Classical bulges in spiral galaxies are known to rotate but the origin of this observed rotational motion is 
not well understood. It has been shown recently that a low-mass classical bulge (C1B) in a barred galaxy can acquire 
rotation from absorbing a significant fraction of the angular momentum emitted by the bar. 

Aims. Our aim here is to investigate whether bars can spin up also more massive CIBs during the secular evolution of 
the bar, and to study the kinematics and dynamics of these CIBs. 

Methods. We use a set of self-consistent N-body simulations to study the interaction of CIBs with a bar that forms 
self-consistently in the disk. We use orbital spectral analysis to investigate the angular momentum gain by the classical 
bulge stars. 

Results. We show that the CIBs gain, on average, about 2 - 6% of the disk’s initial angular momentum within the 
bar region. Most of this angular momentum gain occurs via low-order resonances, particularly 5 : 2 resonant orbits. A 
density wake forms in the C1B which corotates and aligns with the bar at the end of the evolution. The spin-up process 
creates a characteristic linear rotation profile and mild tangential anisotropy in the C1B. The induced rotation is small 
in the centre but significant beyond ~ 2 bulge half mass radii, where it leads to mass-weighted V/cr ~ 0.2, and reaches 
a local Vmax/u"in ~ 0.5 at around the scale of the bar. The resulting V/a is tightly correlated with the ratio of the bulge 
size to the bar size. In all models, a box/peanut bulge forms suggesting that composite bulges may be common. 
Conclusions. Bar-bulge resonant interaction in barred galaxies can provide some spin up of massive CIBs, but the process 
appears to be less efficient than for low-mass CIBs. Further angular momentum transfer due to nuclear bars or gas 
inflow would be required to explain the observed rotation if it is not primordial. 

Key words, galaxies: bulges - galaxies: structure - galaxies: kinematics and dynamics - galaxies: spiral - galaxies: 
evolution 


1. Introduction 


Classical bulges (hereafter CIBs) are the central building 
blocks in many early-type spiral galaxies. CIBs might have 


of cosmic evolution ( 

Kauffmann et al. 

199-4 

1996: Hopkins et aji 

2009; Naab et al. 

2014J) 


; Baugh et all 

or through a 


number of other mech anisms su ch as monolithic collapse of 
primordial gas clouds (lEggen et al.l . Il962l) , the co alescence 
of giant clumps in ga,s-rich primordial galaxi e s (iNoguchil . 
1999; llmmeli et all 12004 lElmegreen et all 200811 . ~vi 
olent disk instability at high-re dshift, dCeverino et al 


vi- 


201511 . multip l e min or mergers (jBournaud et all 12007 : 


Hop kins et al.l. l2010f). and accret ion of small companions 


or satellites ( Aguerri et al.l . 1 200 ill . Although most of these 
studies do not provide quantitative predictions for the 
bulge kinematics, it is generally believed that CIBs formed 
through these processes have lo w rotatio n compared to the 
random motion. For example, iNaab et al.1 (1 20141) showed 
that spheroids produced by minor and major mergers 
(which include CIBs) in full cosmological hydrodynamical 
simulations have a wide range of rotational properties, w ith 
the m assive ones having V/a less than 0.5. lElmegreen et al.1 
(2008) reported dispersion dominated clump-origin CIBs 


with upper limit on V/a ~ 0.4 —0.5, where V is the rotation 
velocity and a is the cen t ral ve locity dispersion. A similar 
study bv llnoue fe Saitohl (120121) suggests that clump-origin 
bulges have exponential like surface density profiles and 
r otate rapidly with V/c t ~ 0.9, resembling pseudobulges 
(iKormendv fe Kennicuttl 120041 ). However, using cosmolog¬ 
ical hydrodynamical simulations with continuous gas accre¬ 
tion, ICeverino et al.l (120151) showed that massive classical- 
like bulges with non-zero angular momenta are produced at 
high redshift but provided no estimate on the bulge V/a. 
Overall, there is a lack of clear quantitative picture of the 
rotational motion induced during the formation of classical 
bulges in numerical simulations. 


Various observational measurements have confirmed 
that CIBs in spiral galaxies possess rotation about their 
minor axis and in most cases in the same sense as the disk 


rotat es ( Kormendv fc Illingworth! Il982t ICannellari et all 
120071 : iFabricius et al.l . l2012h . Disk galaxies both barred 
(e.g., NGC 1023, NGC 3992) and unbarred (e.g., NGC 
4772, NG C 2841) host CIBs with a wide range of masses 


and sizes IKormendv fc Illingworth L 1982^ 
lLaurikainen et all 1 2(1 

cently, Fabricius et al 


Cappellari et al 


KormendvL 119821 


20071 ). More re- 


have obtained detailed kine- 
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matic observations of a large sample of bulges including 
CIBs (whose classification is based primarily on the bulge 
Sersic index). From these studies it is known that CIBs ro¬ 
tate typically with VIa close to an oblate isotropic rotator 
model (jBinnevl . Il978fl . i.e., faster than low-luminosity el¬ 
liptical galaxies, but more slowly than pseudobulges. The 
origin of such rotational motion observed in CIBs remains 
unclear. 


Recent cosmological hydrodynamical simulations which 
include feedback and smooth accretion of cold gas through 
cosmic filaments show that exponential disks could have 
assembled around merger-built C IBs and grown t h rough 
the galaxy’s assembly history (iGovernato et al.1 . l2007t 
lAgertz et all 120111 : iBrook et all l2012h . Once such a disk 
becomes massive and dynamically cold, its self-gravity may 
induce a bar instability leading to rapid formation of a 
bar which would_then int erac t wit h th e preexist ing C1B 
(Hernauis^^Weinberd, Il992t lAthanass oula fe Misiriotis . 
l2002t Saha et all 120121) . In particular. ISaha et al. ( 2012 ) 
showed that an initially non-rotating low-mass C1B could 
absorb a significant fraction of the angular momentum 
emitted by the bar and end up rotating within ~ 2 Gyr. 
However, it is not clear yet how more massive CIBs would 
react to the bar. 

The goal of this paper is to investigate whether angular 
momentum transfer from the bar is an important process 
for the origin of angular momentum in CIBs embedded in 
spiral galaxies. For this, we use collisionless N-body simula¬ 
tions to study the secular evolution in a set of model galax¬ 
ies with initially non-rotating CIBs with different masses 
and sizes. Our simulation results show that, indeed, mas¬ 
sive CIBs could be spun up in their outer parts, typically 
in less than half of the Hubble time. In these models, a 
box/peanut bulge always forms alongside the secular evo¬ 
lution, suggesting that composite bulges may be common. 

The paper is organized as follows. In Section [2] we de¬ 
scribe the initial galaxy models and the set-up for the N- 
body simulations. Section [3] describes the growth of a bar 
and its size evolution in these models. The transfer of angu¬ 
lar momentum and the orbital analysis to account for the 
resonant trapping are described in Section^ The kinematic 
properties of the these bulges are described in Section [5] 
Finally, a discussion of the results and the conclusions from 
this work are presented in Section [6] 


2. Initial galaxy models with massive bulges 

We present here a set of 7 galaxy models c onstructe d using 
the self-consistent method of [Kuiiken fc Dubinskil ( 1995 ). 
Each galaxy model consists of a live disk, dark matter halo 
and bulge. The initial disk is modelled with an exponen¬ 
tially declining surface density in the radial direction with 
a scale-length Rd, mass Md and with a sech 2 distribution of 
stars with vertical scale-height h z . The initial radial veloc¬ 
ity dispersion of the disk stars follows an exponential dis¬ 
tribution with a radial scale length R a = 0.5 Rd such that 
the disk has a constant scale height throughout. The live 
dark matter halo is modelled with a lowered Evans model 
and the classical bulge with a King model. The dark mat¬ 
ter halo has a core in the inner regions of the galaxy model 
and produces a nearly flat rotation curve in the outer parts 
of the disk. For detail s on constructing these models, the 
readers are referred to ISaha et al . (|2010L [2012 b 



R (kpc) L 3?G bulg j 


Fig. 1 . Initial total circular velocity curves (a), Toomre Q 
profiles (c), the total circular velocity curves at 4 Gyr (b) 
and Mass-size relation for the initial CIBs (d). 


We scale each model such that the initial disk scale 
length R d = 4 kpc and the circular velocity at ~ 8.5 kpc is 
~ 200 kms -1 . The initial circular velocity curves for all the 
models are shown in Fig.QJa). Each model has a character¬ 
istic profile for the Toomre Q(R) = cr r (i?)«;(i?)/3.36GE(i?), 
where ay, n and E denote the radial velocity dispersion, 
epicyclic frequency and the surface density of stars; see 
Fig. HTc). The stellar disks in our sample range from dy¬ 
namically cold to fairly hot. Tabic [T] summarizes the initial 
parameters of the disk, bulge and dark halo. 

The bulges in our galaxy models are initially non¬ 
rotating and dispersion dominated. The Sersic indices for 
these King model bulges are generally close to 1 or less. 
Recent findings by show that barred galaxies can have clas- 
sic al bulge like components with sersic ind i ces cl ose to 1 or 
less lLaurikainen et al.l (|2007l) : lErwin et all (|2015l) . The den¬ 
sity profiles in the inner region are flatter than i? 1 / 4 profiles 
which are known to represent the surface brightness profiles 
of traditional CIBs (Ide VaucouleursL 119531 : iFisher fc Drorvl 
I2010D . Unlike i? 1 / 4 profiles, the density profile of a King 
model bulge has a well-defined outer boundary, defined by 
a truncation radius Rb given in table [T| and other well- 
studied properties which motivated us to use them as a 
simple model for the bulges. Setting up of the stellar kine¬ 
matics for the bulge stars is relatively straightforward since 
it is modelled by an analytic DF (King model, fb{E)) that 
depends only on the energy integral (E). Assigning a veloc¬ 
ity to a bulge star is done in the following way: correspond¬ 
ing to a star’s position, we find the local maximum of the 
DF and then employ acceptance-rejection technique to find 
a velocity. The 3 components of a star’s velocity (v x ,v y ,v z ) 
are randomly selected from a velocity sphere with radius 
equal to the local escape velocity. The velocity ellipsoid is 
isotropic by construction. The DF accepts 3 free parame¬ 
ters of which, Oh determines the velocity dispersion of the 
bulge stars; higher the <Jb larger the value of the velocity 
dispersion. Modelled in this way, our initial CIBs are kine¬ 
matically hot spheroidal stellar systems that are similar to 
typical observed CIBs. 
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Table 1. Initial bulge, disk and halo model parameters, 
ordered by the ratio of C1B to disk mass. Column (1): model 
name. (2): Ratio of bulge to disk mass, Mb/Md- (3), (4): 
tidal radius Rb and half-mass radius Rb, 1/2 of the initial 
C1B, normalized by disk scale-length Rd . (5): Q parameter 
(at 2.5 Rd)- (6): disk mass. (7): mass of dark matter halo 
Mh, in units of Md . (8): contribution of bulge and disk 
together to the total circular velocity at 2.2 R d . 


Models .gk 


Rb. 

Rd 


x b,l/2 

Rd 


Q 


M d 

( 10 % 


Mh 

M d 


c.bd 



RCG076 

RCG110 

RHG109 

RCG101 

RCG100 
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Fig. 2. Face-on surface density maps for all model galaxies 
after 4 Gyr of evolution. 


To further compare our initial CIBs with observed CIBs, 
we performed a 2D bulge-disk decomposition of our initial 
galaxy models. We first made FITS files for al l the galaxy 
models and analysed them with GALFIT (|Peng et al.l . 
2002 :. In Fig. EJd) we show the effective radii and masses 
for our initial CIBs, and compare them with the relation 
R P on Mu with a = 0.3 for observed CIBs found bv lGadottil 
(f2009t) . Most of our initial CIBs follow this relation closely, 
while models RHG109 and RCG076 are lcr away from the 
mean relation. 

We evolve each galaxy model in isolation to examine 
the evolution of the bulge shape, morphology and kinemat- 
ics. The simula t ions a re performed using the Gadget code 
(|Soringel et al.l . l200lh which uses a variant of the leapfrog 


= i (Gyr) 

Fig. 3. Time evolution of the bar strength in our galaxy 
models with pre-existing classical bulges. 


method for the time integration. The forces between the 
particles are calculated using the Barn es fc Hut (BH) tree 
with some modification (ISpringel et all 1200111 with a toler¬ 
ance parameter Q to i =0.7. We compute the virial equilib¬ 
rium condition of a galaxy model at different times using 
the Gadget output and found that some of the models were 
although initially out of equilibrium (by about a few per¬ 
cent) settles down close to virial equilibrium (by about a 
percent or less) within a rotational time scale. A total of 
2.2 x 10 6 particles is used to simulate each model galaxy of 
which 1.0 x 10 5 are in the C1B, 1.05 x 10 6 in the disk and 
1.05 x 10 6 in the dark matter halo. The softening lengths for 
the disk, bulge and halo are all unequal and are chosen so 
that the maximum force from parti cl es o f all s pecies (bulge, 
disk, halo) is nearly the same (iMcMillan fc Dehnenl . 120071) . 


3. Growth of a bar and its size evolution 

Although ab out 70% of disk galax i es in the local univers e 
are barred (|Eskridge et all 120001 : iBarazza et al.l . 120081) . 
the formation of a bar is still not fully understood. The 
growth rate of a bar in a disk galaxy depends on var¬ 
ious parameters of the disk and its gravitational inter¬ 
action with surrounding the dark matter halo and pre¬ 
existing C1B (if present). Numerous studies of N-body 
simulations show that a bar forms and grows rapidly in 
a cold rotat i ng stellar disk with Toomre Q close to 1 

(Iff nW 
1200 


1971: Sellwood fc Wilkinsonl. 119931 : lAthanas soulal 


_ Dubinski et al., 2009, and references therein). Swing 

amplification and cooperation of orbital 

streams (iEarn fc Lvnden-Belll ll99fil) are thought to play 
a key role in the bar growth. The linear growth of a bar 
is directly affected by the Toomre Q of the stellar disk as 
can be seen from Fig. [3] A bar grows rapidly in a cool disk , 
whereas it grows rather slowly in a hotter disk (ISaha et all 
l2010f) . The slow growth of a bar occurs through non-linear 
processes; in particular, the resonant gravitational interac¬ 
tion with surrounding dark matter halo is known to play a 
major role. This h as been in vest igate d by s evera l aut hors 
in th e past, dDeb attista fc Sellwo odl. 1998 : Athanassoula. 


2002 


2007 


iHollev-Bockelmann et ah . 20051: IWeinb erg fe Katz . 
ICeverino 


fc Klypin. 20071 ISaha fc NaabT 20131).' 


The role of a pre-existing C1B alone on the formation 
and evolution of a bar has not been fully investigated. Bar 
formation could, in principle, be hampered by a central 
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Fig. 4. The upper panels show the distribution of bulge stars with frequency (17 — 17 b)/k at different times throughout 
the secular evolution in model RHG097. The lower panels show the net change in the angular momentum of the selected 
stars with respect to the previous time. The vertical dotted lines indicate the most important resonances: from left they 
are — 1 : 1, 4 ; 1, 3 : 1, 5 : 2 and 2 : 1. As time progresses, more stars are trapped by the ILR of the bar in the stellar 
disk. However, most of the angular momentum transfer occur through the 5 : 2 resonance. 


mass co ncentrat ion (Hasan et allll993llSellwood fc Moorel 


19991; lAthanassoula et alT 20051b because the bulge could 
cut the feedback loop required for swing amplification by 
placing an ILR (inner Lindblad resona nce) near t he cen - 
ter of a galaxy. In previous studies fe.g.. lSaha et al.l . l2012t) . 
a low mass initially non-rotating C1B facilitated the bar 
growth by absorbing angular m omentum through the ILR. 
Later, ISaha fo Gerhard! (20131 ) investigated the effect of 
varying degree of initial rotation in the preexisting C1B. 
But a detailed investigation of the effect of massive bulges 
on bar formation is yet to be done. 


All the galaxy models in our simulation sample have a 
preexisting C1B of different mass and size (see Table [U, 
and they all form a bar in their stellar disks. Fig. [2] shows 
surface density maps for all stars in each galaxy model at 
the end of 4 Gyr. The models are arranged such that the 
initial value of the bulge-to-disk mass ratio increases grad¬ 
ually from the top left corner (RCG076) of Fig. [2] to the 
bottom right corner (RHG097). Note that Toomre Q alone 
is not the factor deciding the final outcome of these models. 
Dark matter distribution and the bulge might also play an 
important role. Most of these models do not form long-lived 
two-armed spirals. 

In Fig. [3l we show the time evolution of the bar ampli¬ 
tude (A 2 /A 0 ) in all our models. A 2 and Aq are the m = 2 
and m = 0 Fourier component of the surface density respec¬ 
tively. The fastest bar growth is seen in model RCG100 
(Q=1.4). From Fig. Oja) and Fig. [3l we see that gener¬ 
ally the steeper the slope of the initial rotation curve, the 
faster the growth rate of a bar (see model RCG100 versus 
RCG099 and RHG076). However, the correspondence be¬ 
tween growth rate and slope of the rotation curve is not 
one-to-one. For instance, the linear growth rate of the bar 
in RCG076 is higher than RCG110 although the latter has 
a steeply rising rotation curve. Faster bar growth is also 
supported by higher values of (see Table [TJ. Note the 
bar in model RHG109 (Q=2.8) is weakest of all. 


Once formed, the size of a bar increases as it ex¬ 
changes angular momentum with the surroundin g dark 
matter halo and the pre-existing bulge (Saha et all 
120121) . In our simulations, we calculate the bar sizes 
at different ep o chs using the method explained in 
[Saha fe Gerhardl (|2013l ). As these bars grow and become 
self-gravitating, they undergo the well known buckling in- 
stabilit y which leads to th e forma t ion of a boxy/peanut 


bulge 

(Combes & Sandera 1981: Pfenniger & Norman 

199C 

; Raha et al! 1991; Martinez-Valouesta & Shlosman 

2004 

)• 

The galaxy models with prominent boxy/peanut 


bulges end up with a peak in the rotation curve in the boxy 
bulge region, see Fig.[T[b). The rotation velocities shown in 
Fig. [l](b) are calculated at the end of 4 Gyr; by that time 
nearly all the bars are fully grown and have gone through 
the buckling instability. The final bulges in our simulations 
are thus composite bulges which are a superposition of a 
C1B and a boxy/peanut bulge formed from the disk stars. 


4. Angular momentum transfer and orbital analysis 

A bar grows by trapping more and more orbits in its 
main orbital families (e.g., xl orbit family). While it does 
so, the stars loose angular momentum to the dark mat¬ 
ter halo and the C1B in the model. Since the distribu¬ 
tion function, f(E), of the initial C1B is described by a 
King model where it is a function of energy alone, the 
gain of angular momentum by the bulg e at a given res¬ 
onance is always positive (Saha et al! 120121) . The dark 
matter haloes in our model galaxies also gain angular 
momentum but here we are concerned with investigating 
the structural and kinematic changes in the C1B that are 
brought about by the angular momentum gain. To under¬ 
stand the impact of a bar on the pre-existing C1B, we com¬ 
pute the angular momentum transfer between the bar and 
the bulge in our simulat i on using orbital spectral analysis 


feinn e v fc Spergel 1982: Marti noz Valpncsla & ShlosmanL 

2004 iMartinez-Valpuesta et all . 120061) . We focus. in partic¬ 
ular, on two models RHG097 and RCG076 which have the 
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Fig. 5. Same as in Fig. [3] but for model RCG076. 


most massive and the least massive CIBs in our sample, 
with Mb/Md = 0.43 and 0.11 respectively (see Table [0. 

Fig. U shows the results of the orbital spectral analysis 
for the bulge stars in model RHG097 at different times dur- 
ing the evolutio n. As in the case of a low-mass C1B, (see 
ISaha et all 1201 2ll . more and more bulge stars get trapped 
by the 2 : 1 resonance of the rotating potential. However, 
unlike the low-mass case, not much angular momentum is 
gained through the ILR of the bar by this C1B. In fact, at 
~ 2 Gyr (2nd panel in Fig. [4]), the bulge stars are seen to 
be loosing angular momentum through the ILR and this 
continues subsequently with gradually diminishing magni¬ 
tude (except at ~ 4 Gyr, fourth panel from left). In essence, 
we find that the massive GIB in RHG097, on average, loses 
angular momentum through the 2:1 resonance. However, 
and interestingly, we find that a lot of angular momentum 
is gained via 5:2 resonant orbit families (corresponding to 
the peak around (H — f lb)/K = 0.4 in Fig. 0, and this oc¬ 
curs consistently throughout the evolution. In N-body bars, 
these 5 : 2 orbits are known to be stable periodic orbits like 
t he 3 : 1 an d 4 : 1 families inside the corotation of the bar 
(|Voglis et all . l2007lf . 

Fig. [5] shows the results from the spectral analysis of 
model RCG076. We found a similar trend as in model 
RHG097 - a fraction of C1B stars in both models is grad¬ 
ually trapped at the 2 : 1 resonance of the bar. This per¬ 
sists in all CIBs in our simulations (see, Tabled]). However, 
the details and the dominant mode of angular momentum 
transfer vary. At t = 1.96 Gyr (second panel from left of 
Fig. [5], the initial bulge of RCG076 loses angular momen¬ 
tum through the 2 : 1 resonance. But at subsequent times, 
it gains angular momentum consistently through 2 : 1 
orbits. This C1B also gains angular momentum through 
other resonances: during the early phase of bar growth, 
it gains angular momentum through the 3 : 1 resonant 
orbits but gradually, the peak shifts from 3 : 1 to 5 : 2 
resonance (see the right-most panel in Fig. 0. It is in¬ 
teresting to notice that these resonances are rather broad 
compared to the usual 2 : 1 resonance - a feature that is 
common in bo th the bulg e s of R HG097 and RCG076 as 
well as that in ISaha et al. ( 2012 b In other words, associ¬ 
ated with these resonances (e.g., 3 : 1, 5 : 2), there are 
many off-resonant particles that contribute to the net an- 


- 0 . 15 - 0 . 10 - 0.05 0.00 0.05 0.10 0.15 



Fig. 6. The l = 2, m = 2 density wake mode in the C1B of 
RHG097. Dashed lines in all panels indicate the location of 
the bar in the stellar disk while solid lines denote the phase 
angle corresponding to the ’bar’ mode in the classical bulge 
(GIB). Towards the later phases of the evolution, the CIBb 
becomes progressively stronger. The red arrow denotes the 
rotation of the disk bar. Note that at 1 Gyr, a clear bar has 
not developed yet. The unit of time is in Gyr. 


gular momentum transfer. While the two CIBs, in models 
RHG097 and RCG076, gain angular momentum through 
5 : 2 resonance, we do not see any angular momentum 
transfer through the —1:1 orbital families outside coro¬ 
tation whic h contri buted in the case of the low-mass C1B 
studied by ISaha et ail ( 2012 1 using 10 million particles. It 
is not clear to us whether this is due to comparatively 
lower particle resolution in the current models (2.2 million 
particles). Since low resolution noisy simulation could ar¬ 
tificially enhance the sta r-star enco un ters and kno c k star s 
off their resonant orbits, iHollev-Bockelmann et all ( 2005 1: 
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iDubinski et al.l (2009) based on their explicit study on the 
model resolution and resonant behaviour, suggested that 
models with a few million particles generally show conver¬ 
gent behaviour. 

Other than resonant orbit families, there is also a con¬ 
tinuous (chaotic orbits) component of angular momentum 
transfer of comparable importance in both models RHG097 
and RCG076 (around O — D p /k ~ —0.3 in FigsU and [5]). 
We have changed the numerical resolution to calculate the 
potential on a finer grid and repeated the above analy¬ 
sis to rule out numerical artifacts. Also ISaha et ahl (12012H 
found significant contributions to the net angular momen¬ 
tum gain of the C1B via non-resonant chaotic orbits. For 
the two bulges in RHG097 and RCG076, the sign of the 
angular momentum transfer changes during the evolution. 
They both gain angular momentum initially via the non¬ 
resonant orbits, but in the later phases these bulges mostly 
loose angular momentum via the stochastic orbits. Over the 
entire period of evolution, the non-resonant stochastic or¬ 
bits cause a net loss of angular momentum from these bulge. 
However, we have verified that the total angular momentum 
gain through resonant and non-resonant orbits was positive 
at all times. In the following, we look at the density wakes 
created in these bulges by the bar. 

4.1. Density wakes and bulge harmonics 

In this section, we examine the density wakes in the CIBs, 
primarily the l = 2, m = 2 spherical harmonic modes, as the 
bulge stars interact with the bar. Fig. [G] shows the density 
perturbation corresponding to the l = 2, m = 2 spherical 
harmonic mode projected onto the equatorial plane of the 
bulge at different times during the evolution of RHG097. 
Such a density wake arises naturally as the bulge stars 
interact gravitationally with the rotating bar potential. 
If dark matter particles interact with the stellar bar, a 
so ca lled "halo -bar’ gets create d in the dark matter halo 
(iPebattista fc Sellwoodl . l200f)l : iHollev-Bockelmann et al.l . 


- 0 . 15 - 0 . 10 - 0.05 0.00 0.05 0.10 0.15 



Fig. 7. Same as in Fig. [0] but for RCG076. 


20051 Sellwood fc Debattista . 2Of) fit Athanassoulal . 120071 : 
Weinberg fc Katd . 12007 : ISaha fc Naabl l2013lh In the 


present case, we see that the perturbed density in the GIB 
corresponding to this l = 2, m = 2 mode becomes stronger 
as time progresses, in compliance with the increased num¬ 
ber of trapped bulge stars at the 2 : 1 resonance (upper 
panels of Fig. [I] and Fig. [5]). The disk bar rotates anticlock¬ 
wise and its instantaneous position angle is denoted by the 
dashed line (see Fig. ©. 

Initially, the density wake appears in the outer parts 
of the C1B of RHG097 in the form of a spiralish feature 
(at about 1.5 Gyr when the disk bar forms) and become 
stronger eventually taking the shape of a bar-likc mode, 
although not as clear as in the model RCG076 (shown be¬ 
low). The phase of the l = 2,m = 2 mode in the bulge 
of RHG097 reveals that the disk-bar and the bulge-bar are 
nearly aligned with each other during the entire evolution. 
Although there are instances when the two bars are clearly 
misaligned by a very small angle (mostly in the outer parts, 
see Fig. Such_a mi salignment can cause dy namical fric¬ 
tion (IChandrasekhaii 1 19431 iTremaine fe Weinberd . I1984D . 
thereby leading to the transfer of angular momentum ei¬ 
ther from the disk-bar to the bulge or vice-versa. Fig. Q] 
shows that there is little angular momentum exchange via 
the 2 : 1 resonance probably because of the small angle 
misalignment between the two bars of RHG097. 

A similar analysis as above is shown in Fig. [7] for the 
low-mass bulge model in RCG076. The development of a 
bulge-bar through the l = 2,m = 2 mode follows a sim¬ 
ilar evolution as in model RHG097, albeit with some dif¬ 
ferences. As time progresses, the bulge in RCG076 grows 
a comparatively more pronounced bulge-bar. Initially, the 
density wake lags the disk-bar. At t = 1.48 Gyr, a part of 
the bulge-bar appears to be leading the disk-bar. However, 
during the later phases of evolution, this bulge-bar lags be¬ 
hind the disk-bar and gains angular momentum through 
the l = 2, m = 2 mode, in compliance with the angular mo¬ 
mentum gain via 2 : 1 resonant orbits (Fig. [5]). So in this 
case, we have a consistent picture of angular momentum 
transfer via the 2 : 1 resonance and the misalignment of 
bulge-bar and disk-bar. 

Inspection of Fig.[G]and Fig.Qshows that the peak of the 
l = 2, to = 2 density wake develops at different spatial loca¬ 
tions in the two bulges - for the massive GIB (R.HG097) it is 
pronounced in the outer parts while for the low-mass case 
(RCG076) it is confined to an intermediate radial range. 
Thus, only the stars in the outer parts of the GIB in model 
RHG097 take part in the rotational motion; see the line- 
of-sight velocity maps of Fig. [8] It is likely that both figure 
rotation and streaming motion within the bulge-bar con¬ 
tribute to the net rotational motion. 


5. Bulge kinematics 

The orbital configuration of the CIBs in our model galax¬ 
ies changes as a result of the angular momentum gain via 
the bar-bulge interaction, and this must manifest itself in 
the kinematic properties of the CIBs. In this section, we 
focus on the stellar kinematics, in particular on the rota¬ 
tional motion of the classical bulge component induced by 
the angular momentum gain. By separating the GIB parti¬ 
cles from the boxy/peanut bulge that forms from the disk 
through the bar buckling instability, we investigate the de- 
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X/R d X/R d 

Fig. 8. Surface density (upper panel), line-of-sight veloc¬ 
ity and velocity dispersion maps (2nd to 4th panel) of the 
C1B in RHG097 at different times during the evolution. 
These images are taken at 90° projection (edge-on view) 
and the major axis of the bar is aligned with the X-axis. 
Clear signatures of rotation are seen at 4 Gyr. Color bar at 
top represents density, middle the velocity and bottom the 
velocity dispersion. 


pendence of the acquired rotational motion in the CIBs on 
various parameters. 

5.1. Rotation and dispersion profiles 

In Fig. [5J we show the evolution of the massive C1B model 
RHG097 through surface density, line-of-sight velocity and 
dispersion maps in edge-on projection. As the bulge gains 
angular momentum, it starts picking up rotation. Signs of 
systematic rotation are obvious at t = 2 Gyr, and at t = 
4 Gyr, the C1B rotates faster, with (V/cr) ~ 0.5 in the 
outer parts. It is interesting to notice that the inner parts 
(R < 2R b i/ 2 ) of this C1B do not show much rotational 
motion - this mostly confined to the outer parts of the C1B. 
At the end of 4 Gyr, the initially non-rotating C1B has also 
become slightly rounder and hotter at the center. 

In Fig. [9l we show radial profiles of rotational velocity, 
V(R), velocity dispersion, a(R), and local V(R)/o- m , where 
<Ti n denotes the average velocity dispersion in the central 
region (within 0.5/4, 1 / 2 ), for four CIBs after 4 Gyr of evo¬ 
lution. We illustrate models RHG097, RCG100, RCG101 in 
order of decreasing bulge-to-disk mass ratio, and the rela¬ 
tively compact bulge model RHG109. The slope of the rota¬ 
tional velocity is approximately the same for all four bulges 
shown in the figure. This is also true for the remaining CIBs 
in our sample. The mean velocity dispersion profiles are also 
nearly identical for all these bulges except RHG109 which 
has lower velocity dispersion throughout. The V/<j- m pro¬ 
files for the four bulges are shown at the bottom panel of 
Fig. [9] In all cases, rotational motion becomes significant 
in the outer parts of the bulge, for R > 0.5 Rd — 2Rb, 1 / 2 - 

In some of the models in which the bar grows rapidly, 
reaching its peak amplitude (A 2 , ma x) within a billion years 
(e.g., RCG100, RCG101), the bulge stars acquire significant 


rotation already within 2 Gyr. Thereafter, over the next 

3 — 4 Gyr, the rotational profiles in these bulges change only 
insignificantly. On the other hand, in models that grow their 
bar rather slowly, reaching peak amplitude over 1 — 3 Gyr, 
the rotation velocity in the CIBs continues to change until 
the bar growth approximately saturates at around 4 Gyr. 

Comparing the rotational profiles of our simulated CIBs 
with those in observed barred galaxies is not straightfor¬ 
ward because the kinematic data in external galaxies in¬ 
clude both the C1B stars and the stars in th e disk and a 
possible boxy bulge. Recently, iFabricius et al. ( 2012 1 anal¬ 
ysed a large sample of galaxies containing CIBs residing in 
both barred and non-barred galaxies. The rotation profiles 
of the GIB host galaxies in their sample have a wide range 
from shallow rising profiles, to steeply rising profiles often 
associated with steeply falling velocity dispersion. Visual 
comparison of the V/cr profiles of their bulges with our sim¬ 
ulated CIBs suggest that some of the shallow-rising profiles 
could be similar to our spun-up CIBs. The amplitude of 
the rotational profiles ( V/cr{R )) for the GIB host galaxies 
re ached within their bulg e radius is ~ 0.5 — 1.0; see Fig. 13 
of lFabricius et al. (201.2!). But note that their bulge radius 
is not simply related to the effective radius R e ; for their 
entire sample, it is typically ~ 2 R e . The V/o{R) reached 
by the rotation profiles of our spun-up CIBs are ~ 0.2 at 
2/4,i/ 2 - Beyond this radius, V/o(R) for our CIBs keeps in¬ 
creasing, reaching V/a{R) ~ 0.5 at ~ 4Z4 ,i/2- 

However, the final bulge in our models is always a com¬ 
posite bulge composed of the preexisting C1B and the boxy 
bulge formed from the disk. Fig. [10] shows the line-of- 
sight velocity, dispersion and surface density profiles of the 
composite bulge (C1B and boxy bulge combined) of model 
RHG097. The inner parts of this galaxy model show a mod¬ 
erately rising rotation profile and declining dispersion. The 
size of the boxy bulge is about 0.6/4>ar — 0.6 Rd, not much 
greater than ~ 274,1/2 of the GIB, which is ~ 0.5/4- At 
this radius, the local V/o ~ 1.0. The slope of V/a(R) pro¬ 
file for this co mposite bulge resemble s many of the observed 
pseudobulges (IFabricius et al.l . l2012f l. A more detailed com¬ 
parison with observations is outside the scope of this paper. 

Finally, we report here about the three-dimensional na¬ 
ture of the induced rotational motion in our CIBs. None of 
the massive bulges in our sample rotates cylindrically after 

4 — 5 Gy r, contr ary t o the c ase of the low-mass GIB investi¬ 
gated bv ISaha et al.l ( 2012 1. Evidently, the orbital changes 
that occurred in the high-mass bulges as a result of their 
interaction with the bar are weaker than in the low-mass 
case. 

5.2. Induced rotational kinetic energy and anisotropy 

In Fig. O we noticed that V/cr rises in the outer parts of 
the CIBs whereas the central parts of the bulge remain with 
small rotation (V/cr < 0.1). This can also be readily seen 
from the velocity maps of the C1B of RHG097 (see Fig. [5]). 
Both these facts have prompted us to investigate the actual 
rotational energy that is being imparted by the bar to these 
bulges. The random kinetic energy of a group of stars mov¬ 
ing in the equatorial plane of the bulge, using cylindrical 
coordinates, is given by 

T’rand(-Rj) = T Ian ^ R (Rj) + T Iand ^(Rj), (1) 

where, 
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Fig. 9. Rotation, velocity dispersion, and local V/a radial 
profiles for the four CIBs in models RHG097, RCG100, 
RCG101 and RHG109. The blue and black arrows show the 
half-mass radii for the two bulges of RHG109 and RHG097, 
respectively. 


-i 

T r and,ii(-Rj) = o y ]rnj[vR{i) -v R ] 2 , (2) 

Z i=l 

and 

1 Nj 

Trand ^ ' 'TYlil'Vcpii') Vip] (3) 

Z i=1 

where Nj is the number of bulge stars in the j th radial 
bin and Rj denotes its radius; rrii is the mass of each bulge 
star, vr and v^ are the mean radial and azimuthal velocity 
of the stars. 

Similarly the rotational kinetic energy is given by 

N i 

T I0 t{Rj) = — rriiVy, (4) 

z j=i 

so that the total kinetic energy of a group of stars moving 
in the equatorial plane is 

(5) 



0 0.5 1 1.5 2 

R/ R d 

Fig. 10. Surface density (top), mean line-of-sight velocity, 
and dispersion profiles (middle panel) for the composite 
bulge of model RHG097. The local V/a profile is shown in 
the bottom panel. 



Fig. 11. Fractional kinetic energy associated with the ro¬ 
tational motion for the four bulges in RHG097, RCG100, 
RCG101 and RHG109 in radial bins. The blue and black 
arrows show the half-mass radii for two bulges RHG109 and 
RHG097 respectively. 

Fig. HD shows the radial distribution of the fractional 
rotational kinetic energy in the final C1B bulges (used 
Eqs. 1 — 5). The figure clearly shows that the inner regions 
(R < 2R bl / 2 ) of these CIBs do not acquire much rotation 
from the bar, with final rotational kinetic energy amount¬ 
ing to about 5% of the total planar kinetic energy after the 
evolution. It is the outer regions which acquire most of the 
angular momentum transferred, where for R > 274,1/2 tyP - 


Ttot(Rj) — T rand (Rj) + T Iot (Rj). 
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Fig. 12. Radial variation of the radial, tangential, and ver¬ 
tical velocity dispersions in the four bulges in RHG097, 
RCG100, RCG101 and RHG109 after 4 Gyr evolution. The 
blue and black arrows show the half-mass radii for two 
bulges RHG109 and RHG097, respectively. 



Fig. 13. V/a versus ellipticity for all the CIBs in our galaxy 
models at the end of 4 Gyr. Initially all simulated CIBs were 
non-rotating. Over-plotted are the V/a values for CIBs from 
Fabricius et al. (2012). 

ically 30 — 40% of the planar kinetic energy is in rotation 
for most of the CIBs. Note a possible explanation for this is 
given in section im In terms of mass, this rotating region 
broadly contains about 40% or less of the total bulge mass. 

In Fig. [T2j we show the radial variation of the radial, 
azimuthal and vertical velocity dispersion (tjR,t t v and a z ) 
for four CIBs at the end of 4 Gyr. Initially, all the bulges 
were isotropic by construction. Of the four CIBs, RHG109 
was spun-up significantly less than the other three after 
4 Gyr of evolution. In all four CIBs, the final azimuthal ve¬ 
locity dispersion is slightly higher compared to the radial 
and vertical dispersions, with the largest effect in the out¬ 
skirts where rotation velocity is also highest. Thus Fig. IT2l 
shows that these spun-up bulges tend to become tangen¬ 
tially anisotropic. 

5.3. V/a for the spun-up bulges 

Here we plot our sp un-up jCIBs on the {V/a, e) plot 
dKormendv fc Illingworthl . 4982 ) to illustrate their rota¬ 
tional properties further. Note that while this is easily done 
in simulations because of the unique ID attached to each 


particle, comparing with observed CIBs is non-trivial, be¬ 
cause it is difficult to separate classical bulge stars from 
stars in the boxy bulge in observations. First, we calculate 
V/a values in the traditional sense. We record the peak ve¬ 
locity, V m , which is taken to be a simple average of the outer 
velocities in the radial range ~ 0.9 — 1.2 Rd (excluding out¬ 
liers such as the point with ~ lOOkrns -1 in model RCG100, 
see Fig. ®. This value is divided by the inner velocity dis¬ 
persion, er in , computed by taking an average within the 
half-mass radius of the C1B after 4 Gyr. The resulting V/a 
values are listed as V max /a- ln in Table [2] 

In addition, we also calculate V/a from the edge-on 
kinematic maps of our evolved CIBs, following iBinnevl 
(2005), and denote these by {V/a) < 2 d> ■ Thus, we com¬ 
pute the surface density-weighted velocity and velocity dis¬ 
persion both in the inner {R < 2R bl / 2 ) and outer parts 
{R > 2R bl / 2 ). Consistent with the results of the last sec¬ 
tion, the {V/a )< 2 d> in the inner parts are small (< 0.1). 
In the rest of the paper, we show the {V/a) <2 d> values for 
the outer parts only. 

We calculate the ellipticities (e) of the CIBs by the 
method of diagonalizing the moment-of-inertia tensor. The 
intrinsic ellipticity of the bulge is then e = 1 — c/a, where 
c/a is the ratio of the minor to major axis. In edge-on pro¬ 
jection (where the inclination angle is 90° w.r.t. the LOS), 
the apparent ellipticity would be equal to the intrinsic el- 
lipticity derived using the moment-of-inertia method. We 
compared these axis ratios with those obtained using IRAF 
ellipse fitting, and the agreement between both methods is 
fairly good. 

Fig. Il3l shows the well-known V/a — e plot for all the 
CIBs in our model galaxies after 4 Gyr. Ellipticities and 
V/a values at this time are also tabulated in Table [2] 4 
of the initially non-rotating CIBs from our simulated sam¬ 
ple have final V max /a ln values which formally put them on 
the oblate isotropic rotator line (iBinnevl . Il978l ). while the 
remaining 3 bulges did not spin up as much. For all our 
CIBs, the outer (V/a) <2 d> values are lower compared to 
the V mSLX /a in . This is because the induced rotation profiles 
in these CIBs have a shallow rise and reach their maxima 
in the outermost parts which contain little mass. The final 
{V/a )< 2 d> may be slightly on the low side in our models, 
however, because the initial CIBs were modelled as King 
profiles in which the density of stars drops to zero rapidly 
in the outer parts marked by the truncation radii and con¬ 
tain fewer stars. In contrast , the observ ed CIBs which fol¬ 
low Sersic profiles r 1 n te.g.- ICaon et al.Lll993l) . with n ~ 4 
have comparatively less steeply falling density and contains 
larger number of stars at the outer parts. Our analysis re- 


Table 2. Bar size and bulge kinematics at the end of 4 Gyr. 


Models 

-^bar 

Rd 

Jh 

^bar 

e 

Knax/ ^in 

{V/a)< 2d > 

{R > 2 R b ,i/ 2 ) 

RCG076 

0.61 

3.42 

0.37 

0.48 

0.29 

RCG110 

0.83 

0.97 

0.33 

0.41 

0.25 

RHG109 

0.90 

0.72 

0.44 

0.26 

0.19 

RCG101 

1.10 

0.96 

0.31 

0.68 

0.26 

RCG100 

1.03 

1.04 

0.29 

0.60 

0.23 

RHG099 

0.94 

1.15 

0.19 

0.42 

0.27 

RHG097 

0.93 

1.31 

0.17 

0.48 

0.30 

RCG004 

0.85 

1.52 

0.28 

1.40 

0.77 
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Fig. 14. Dependence of V/a on the ratio of bulge size to 
bar size of the embedded CIBs at the end of 4 Gyr. Initially, 
all the bulges were non-rotating. The star signs indicate the 
C1B of RCG004, as in Fig. M 


veals that it is the outer part of the massive bulges which 
gain most of the angular momentum transferred by the bar. 
Having these outer parts well populated by stars on stable 
orbits would help resonant transfer of angular momentum 
more effectively or vice-versa. We, therefore, expect angu¬ 
lar momentum transfer and the bulge spin up to be greater 
for less steep bulge profiles as might be the case in observed 
CIBs. 

The two star symbols show the equivalent V/a values 
for the low mass C1B model of RCG004 from ISaha et all 
(2012:). They are considerably higher than for the massive 
bulge models reported here. The C1B in this model had sig¬ 
nificantly lower mass (0.067 Md) but similar Rb,\ii/Rd — 
0.21. Clarifying the origin of this difference will require 
a more extensive parameter study of the evolution of low 
mass bulges in bar unstable disks. 

In addition to our simulated CIBs, we have overplotted 
in Fig. Q31 6 CIBs in barred galaxies (NGC1023, NGC2859, 
NGC2880, NGC3521, NGC3992 and NGC4260) and 7 CIBs 
in non-barred galaxies (NGC4698, NGC4772, NGC3898, 
NGC3245, NG C3031, NGC2775, NGC2841) taken from 
iFabricius et al.1 (|2012t) . In most cases, the CIBs in barred 
galaxies rotate faster than their counterparts in non-barred 
galaxies. Whether in these galaxies the bar might be re¬ 
sponsible, in part, for spinning up these CIBs will require a 
separate investigation. 

Fie, fill shows the dependence of normalized {V/a)* on 
the ratio of bulge size to bar size ( Rb/Rbai )■ For the nor¬ 
malization, we used the relation V/a ~ -\/e/(l — e). For all 
bulges (except RCG076 ) the f inal (V/a )* at 4 Gyr, as well 
as {V/a )< 2 d> following lBinnevi (l2005h correlate tightly with 
increasing ratio Rb/Rb&r ■ This suggests that the dynamics 
of the C1B spin-up depends on the relative size of the bar. It 
is possible that there is an optimal size for the rotating bar 
to provide the most angular momentum transfer to the C1B 
through resonant angular momentum transfer and dynam¬ 
ical friction on the bulge stars (see Section HUD; however, 
this needs additional models with large CIBs. Certainly, 
from the models studied here, the spin-up of CIBs is not 
efficient when the C1B is too compact. 


6. Conclusions and discussion 

We have used self-consistent N-body simulations of disk 
galaxies to study the interaction of massive CIBs with the 
bar that forms in the disk. We found that significant an¬ 
gular momentum is transferred by resonant gravitational 
interaction from the bar to the C1B in all cases, in a time- 
scale of a few Gyr. We have analyzed the angular momen¬ 
tum transfer using orbital frequency analysis and examin¬ 
ing the barred density wake in the bulge. Because in our 
models most of the angular momentum is transferred to the 
outer parts of the C1B, it is likely that the induced rotation 
in Sersic bulges would be higher than in those studied here. 
Our main conclusions from this work are as follows: 

1. Massive classical bulges gain as much specific angular 
momentum through spin-up by the bar as do lower-mass 
bulges. 

2. Most of the angular momentum transfer occurs through 
low-order resonances. In particular, we found that a lot of 
angular momentum can be transferred via 5 : 2 resonant 
orbits which are common in the models studied here. 

3. The bar generates a density wake in the bulge due to 
trapped 2:1 orbits. For lower mass bulges in our sample, 
such a wake is found to be misaligned with the bar and 
causes angular momentum transfer while for higher mass 
bulges the misalignment is not substantial. At later phases 
of evolution, the bulge density wake aligns with the bar and 
evolves into the bulge-bar. 

4. The spin-up process creates a characteristic linear rota¬ 
tion profile such that significant rotation is induced beyond 
~ 2i? M/2 . Typical mass-weighted V/a beyond ~ 2 Rb,i /2 
are ~ 0.2, and the local V miix /a\ n reached at the largest 
radii are ~ 0.5. 

5. Contrary to t h e case of the low - mass CIBs studied in 
ISaha et aLl(l2012H : ISaha fc Gerhardl (1 201 3li . the rotation in¬ 
duced in the C1B is not cylindrical. 

6. The spin-up process also creates mild tangential 
anisotropy in the outer regions of the CIBs studied here. 

7. In all models a box/peanut bulge also forms through 
the bar-buckling instability. The final systems therefore 
have a composite bulge with the C1B superposed on the 
slightly larger box/peanut bulge. This suggests that com¬ 
posite bulges may be common in galaxies. 

The main result of this paper is that also comparatively 
massive bulges are impacted by the angular momentum 
transfer mechanism from the bar, which was previously 
i nvestigated only for low-mass CIBs, both non-rotatin g 
(ISaha et al.l . 1 20 12li and rotating (|Saha fe Gerhardl 120131) . 
In our sequence of models, the spin-up is most efficient 
when the size of the bar is somewhat larger than the C1B, 
but not too large. This suggests that larger effects might be 
expected in the presence of gas, which would favour higher 
pattern speeds and mass inflow. If nuclear bars formed as 
a result of central mass build-up, they could lead to ad¬ 
ditional angular momentum transfer. These effects merit 
future investigation. 
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